Casimir effect of electromagnetic field in iD-dimensional spherically symmetric cavities 
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Eigenmodes of electromagnetic field with perfectly conducting or infinitely permeable conditions 
on the boundary of a _D-dimensional spherically symmetric cavity is derived explicitly. It is shown 
that there are (D — 2) polarizations for TE modes and one polarization for TM modes, giving rise to a 
total of (D — 1) polarizations. In case of a D-dimensional ball, the eigenfrequencies of electromagnetic 
field with perfectly conducting boundary condition coincides with the eigenfrequencies of gauge one- 
forms with relative boundary condition; whereas the eigenfrequencies of electromagnetic field with 
infinitely permeable boundary condition coincides with the eigenfrequencies of gauge one-forms with 
absolute boundary condition. Casimir energy for a D-dimensional spherical shell configuration is 
, computed using both cut-off regularization and zeta regularization. For a double spherical shell 

■ configuration, it is shown that the Casimir energy can be written as a sum of the single spherical 

shell contributions and an interacting term, and the latter is free of divergence. The interacting 
, term always gives rise to an attractive force between the two spherical shells. Its leading term is the 

Casimir force acting between two parallel plates of the same area, as expected by proximity force 
approximation. 
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(]^ ■ I. INTRODUCTION 

This decade has observed a tremendous progress in the research in Casimir effect. One of the motivations comes 
—^ from the wide apphcations of Casimir effect in other areas of physics Besides, successful experimental verifications 
^ of this quantum effect by various groups of researchers have greatly stimulated the research in this area. Before the 
turn of this century, the research in Casimir effect focused on simple geometric configurations such as parallel plates, 
rectangular, spherical and cylindrical geometries, where specific coordinate systems are available for explicit solutions 
of the corresponding energy spectrums. For more complicated configurations such as sphere-plate configuration, 
approximation methods like proximity force approximation have to be employed to approximate the Casimir force. In 
\ this decade, more attention was paid to numerically calculate the Casimir energy or Casimir force between arbitrary 
compact objects. Various methods have been proposed such as the semiclassical method, the worldline method and 
the multiple scattering approach. 
• • , Most of the works in Casimir effect, especially those involve real materials considered (3-1- l)-dimensional Minskowski 
' spacetime. Nevertheless, the exploration of Casimir effect in higher dimensional spacetime has a long history. One 
, of the pioneering works is the work by Ambj0rn and Wolfram [3], where they derived Casimir energy for massless 
5_i scalar field and massless vector field (electromagnetic field) in a Z?-dimensional rectangular cavity. Despite of the 
5^ , long history, the research on Casimir effect in higher dimensional spacetime has been focused on massless scalar field. 
Relatively few works have considered massless vector field since the eigenmodes are harder to solve compared to the 
scalar case. Very often the results for vector field were incorrectly claimed as a simple multiple of the results for scalar 
field. Except for the configuration of a pair of parallel hyperplanes, this claim is not true. In case of D-dimcnsional 
rectangular cavity, this can be seen from the results of . 

Motivated by Casimir 's intriguing model for electron [3], Boyer Q calculated the Casimir force acting on a perfectly 
conducting three dimensional spherical shell and found that it is repulsive. This result has later been confirmed by 
several groups of researchers [5Hl2l| using different methods. In this case, the eigenmodes can be divided into transverse 
electric (TE) modes and transverse magnetic (TM) modes. Except for the absence of the I = modes, the TE modes 
is the same as the eigenmodes for massless scalar field with Dirichlet boundary conditions, and the TM modes 
correspond to the eigenmodes of massless scalar field with Robin (mixed) boundary conditions with a suitable Robin 
parameter. For massive scalar field confined in the interior and exterior of a three dimensional spherical surface, 
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the Casimir energies are computed in [13|. The extension to higher dimensions was first carried out in [ij], where 
the Casimir stress acting on a D-dimensional spherical shell was computed for massless scalar field using Green's 
function method. This was later reconsidered in |15l| using zeta function method. An attempt to generalize this result 
to massless vector field was carried out in fio', '17*]. By analogy with the D — 3 case, the Casimir stress acting on 
a D-dimensional perfectly conducting spherical shell is divided into contribution from TE modes and TM modes, 
where the contribution by TE modes is equated to the contribution by massless scalar field with Dirichlet boundary 
conditions after the omission of the I = mode, and the contribution by TM modes is equated to the contribution 
by massless scalar field with Robin boundary conditions with a suitable Robin parameter. Later in p^ . a systematic 
approach for scalar field, spinor field and electromagnetic field in Z)-dimensional spherical cavity was given, and the 
numerical results for Casimir energies when 2 < D < 9 were tabulated for these different quantum fields. Although 
the results for massless scalar field in D-dimensional spherical cavity is authentic, there are some reasonable doubts 
that can be raised on the results for electromagnetic field. In the work |2|, it was found that the electromagnetic 
field in a Z?-dimensional cavity has (D — 1) polarizations. However in |lq-[l8|. the electromagnetic field only has two 
polarizations in a Z?-dimensional spherical cavitiy. In our opinion, this inconsistency is because the electromagnetic 
field in Z?-dimensional spherical cavity was considered in [16l - ll8l | as a pure direct generalization of the D = 3 case, 
without reference to the equation of motion satisfied by electromagnetic field, as was considered in Q . One of the 
aims of this article is to resolve this problem. 

In a, {D + l)-dimensional Minskowski spacetime with metric 

the field strength of electromagnetic field is represented by the two form F — dA = F/^^^dx^dx'^ . As in we introduce 
the potential 

A = A^dx'' |:=^A^dxM 

\ Ai=0 / 



SO that 



F = dA=(^-^] dx^^dxr 
\dx'' dx" 



The action is the standard one: 



where F*^^ — g'^'^g'^'^F^^, and the equation of motion is 

1 d ( /r^^,.„\ 1 d 



dx 



77 {VHF^" 



\g\ dx- 



dA^ dAn 
dx^^ dx^ 



0. 



(1) 



As is well-known, there is a gauge degree of freedom given by 

Ai^ A + dLp, 

where Lp is any function. This degree of freedom can be eliminated by imposing the radiation gauge 

1 d 



Ao = 0, 



0, 



A^' = gf'-A^ 



(2) 



Under this gauge condition, the equation of motion ([T]) is equivalent to the Laplace equation 



AA„ 



9m 



d 



dx"" 



OA,. dAr, 



dx''i dx'^ 




= 



(3) 



on the one-form A. There has been a work [l9| which computed Casimir energy for Laplace operator acting on 
p-forms (0 < p < D) on a _D-dimensional ball . The boundary conditions imposed on the boundary of the ball - 
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the {D — l)-sphere, are the so-called absolute and relative boundary conditions [20j. In spherical coordinates: 



-r cos 01 

:r sin^i cos 62 



r > 0, 

<e',; < TT, l<i<D 
<</> = Od-1 < 27r, 



(4) 



„D-1 



-r smWi 
-r sin 01 



. sm VD-2 cos ( 
. sin 9u-2 sin (J 



the absolute boundary conditions for one- forms on a D-dimensional ball read as [l9l [21j: 

dAe 



dr 



0: 



0; 



and the relative boundary conditions read as [l9|, l2]J : 



0, 



d_ 

dr 



D - 1 



A, 



0. 



(5) 



(6) 



However, as a potential for the massless vector field, the boundary conditions imposed on A as discussed in are 
completely different. In fact, the boundary conditions are imposed on the field F = dA rather than the potential A. 
Generalizing the perfectly conducting and infinitely permeable boundary conditions in _D = 3 dimension, for any Z?, 
the perfectly conducting boundary condition is specified by [2J: 



^lUl...VO-2 



0, 



and the infinitely permeable condition is specified by 

Here is the spacelike vector normal to the surface S, and 

is the {D — l)-form dual to F. These two boundary conditions are gauge invariant since d{dtf) = 0. In terms of the 
spherical coordinates, the perfectly conducting boundary condition on the boundary of a £)-dimensional ball can 
be explicitly written as 



fdAe^ dAg^ 



36, 



0, 



l<i<i <D~l, 



and the infinitely permeable boundary condition can be written as 

'dAn, dA,. 



dr 80, 



0, 



l<i< D-l. 



(7) 



(8) 



At the first sight, the absolute and relative boundary conditions ^ and ^ and the perfectly conducting and infinitely 
permeable boundary conditions ([7]) and ^ are completely different. A result we are going to show in this article 
is that on a Z?-dimensional ball, the absolute boundary condition is equivalent to the infinitely permeable boundary 
condition, and the relative boundary condition is equivalent to the perfectly conducting boundary condition. This 
gives a physical interpretation to the absolute and relative boundary conditions considered in [19]. 

After writing out explicitly the eigenmodes of the electromagnetic field in spherical coordinates, we compute the 
Casimir energy using cut-off regularization and zeta function regularization. The Casimir energy for the region inside 
the spherical shell and the region outside the spherical shell will be computed separately. The Casimir energy for a 
single spherical shell configuration is obtained as a sum of the two. In the last ten years, there has been an interest in 
considering double spherical shell configuration. In [2^ , the stress energy tensors for a massless scalar field inside and 
outside a Z)-dimensional spherical shell, as well as inside the region between two spherical shells have been computed. 
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This has been generalized to spinor fields and spacetimes with different geometries [23l - l25| . In [26l - l28j . the Casimir 
effect due to electromagnetic field with perfectly conducting boundary conditions on two concentric three dimensional 
spherical shells was considered. In this article, we are going to generalize this to D-dimensions. 

The layout of this article is as follows. In Section |lll we derive explicitly the eigenmodes of electromagnetic field 
in D-dimensional spherically symmetric cavity using spherical coordinates. Both perfectly conducting and infinitely 
permeable boundary conditions are considered. In Section IIIH we compute the Casimir energy inside a spherical 
cavity by cut-off and zeta regularizations. In Section IIVI we show that the Casimir energy inside an annular region 
bounded by two spheres can be written as a sum of three terms: the Casimir energy inside the sphere with larger 
radius, the term that can be interpreted as the (renormalized) Casimir energy outside the sphere with smaller radius, 
and the interacting term. In Section |Vl we consider the (renormalized) Casimir energy outside a spherical shell. In 
Section rvTl we consider the Casimir energy of a single spherical system which is the sum of the Casimir energy inside 
the sphere and the Casimir energy outside the sphere. In Section IVIIl we consider a double spherical shell system 
whose Casimir energy is obtained by summing the Casimir energy inside the sphere with smaller radius, the Casimir 
energy in the annular region between the two spheres, and the Casimir energy outside the sphere with larger radius. 
Using the result of Section IIV[ this can be rewritten as the sum of the Casimir energies of the two spheres, and the 
interacting term. In Section [Villi we present a method to compute the asymptotic behavior of the Casimir force when 
the separation between the spheres is much smaller than the radii of the spheres. A brief conclusion of the results 
obtained in this article is given in Section |lXl 

Throughout this article, we assume that D > 3 and we use the units where h = c — 1. 



II. EIGENMODES OF THE ELECTROMAGNETIC FIELD IN SPHERICAL COORDINATES 



In this section, we will write down explicitly the eigenmodes of the electromagnetic fields in (D + l)-dimensional 
spherical coordinates that satisfy the radiation gauge ([2]). In spherical coordinates (|4]), the metric tensor g^j,^ is 



ds^ — g^^dx^dx'^ = dt^ — dr^ — 



d0f 



(9) 



and 



D-2 



We want to find explicitly the one forms 

A = Ai_,dx'^ = Atdt + Ardr + Ae^dOi 



Ag^_.^d9D-2 + A^dcj) 



which are solutions of the equation of motion ([T} subject to the gauge condition In spherical coordinates, the 
gauge condition is given by 



At =0, 







1 d 



'dxf^ 



{^/\9\A^) 



1 d 

„D-1 ^ 



D-1 



r^-'Ar 



i=l 



sin^ 6-, 



0,Aa 



(10) 



In spherical coordinates, the equation of motion ([1} or the Laplace equation ([3]) equivalent to ([TJ under the radiation 
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gauge pO|) are quite complicated. We will consider another set of equivalent equations. Adding -25 to the Laplace 

r 

equation ([3]) when fi = r, we obtain the equation 



D-l 



d^Ar D + ldAr D-l , sr 



1 



d 



Qj.2 



dr 



^ r2 sin^ 01... sin^ dj-i sin^^^^^ 9^ d9j 



sm 



D-j-l 



OAr 

'99, 



= 0. (11) 
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2 cos Oi 

For 1 < j < -D — 1, adding 7—25 to the Laplace equation (El) when n = 9i, we obtain the equation 

sm Oi 



sm 



' 5/15 fsin--^-0,A. 



r2 sin2 01 . . . sin^ 0,_i sin^ 0^ V sin^-*-^ 9, dO, 



E 



sm' . . . sm' tlj-i sm" ■' " t/j oOj \ oOj 



sm 



W (12) 



2 9^^ „i-4cos0i 1 dAg 



2E 



^ sin 0j r2 sin 0i . . . sin 0j_i d9i 



2COS0. 1 9 . „_i 2cos0, ;-A 1 ^ fsin^-J-i A« ^ - 

sin0, r^-iar ^''^ sin0, sin^ 0^ . . . sin^ 0,_i sin^-^-i 0, 90^ ^^.^ej - U. 

Notice that (fTTj) only contains A^, and (fT2|) only contains Ar and with 1 < j < i. To further simplify the problem, 
we divide the energy eigenmodes of the field F — dA into TE modes and TM modes as in the D = 3 case. Since the 
component Ei of the electric field is identified with Foi and the component Bi of the magnetic field is identified with 
{Fjk ■ ^^j<k<D,j,ky^i} [2], on a ZJ-dimensional ball or annular region, TE modes and TM modes are defined 
as modes satisfying 

TE modes : Fq,. = 0, TM modes: Fg^e^ =0, I < j < k < D -1. (13) 

Interestingly, we find that this division into TE and TM modes are exactly the same as the decomposition 

A = A^ + 2t((^) 
for a one form A satisfying the transversality condition 

1 d 



{■Ma'') 







discussed in [1^, [20]- To be more precise, the condition for TE modes Fp^ = is satisfied if and only if A,. — 0, and 
this is exactly the characterization of in [l^ [l^] ■ For the term 2l(iy9) , notice that the condition on TM modes 
implies that the corresponding A can be written as 

A = Ardr + e~''^*/(r)i, 

where A is a closed one form (i.e. dA = 0) on the unit [D — l)-dimensional sphere 5"^^^. Since Z? > 3, the first 
cohomology of 5*^"^ vanishes. Therefore there exists a function (p on 5*^^^ such that 

d(p ^ A. 

In other words, for 1 < i < — 1, 

The equation of motion (IT|) for /i = 0^ , 1 < i < Z) — 1 becomes 

3s^e^ 5 / rrift^g^ttTj. \ 9e,e, d / rrig.g^rrrp \ _d'^Ag^ 1 d 



This implies that 

Ar = e-'"*5M^, 



j,D 3 Qj, 



r 



3 


'dAr 


dAg^ 






[89, 


dr 


) 
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where 



^ dr ^ dr \ dr 

Since the Laplace operator on scalar functions on the {D — l)-dimensional unit sphere S^^^ is 

D-l 



- - V 1 ^_ A [ sin^— 1 e — 



[n;ilBi 



sm^ ( 



the equation for Ar dTTI) and the gauge condition imply that 



2.^~ ^ + 1 dgjr) , iJ - 1 , , - , gjr) ^ (o) , ^ 



3 f^j. 

These are satisfied if and only if 



d^gir) D + ldg{r) \^ - {D - 1) . . ^ 

-jV + ^ ^ffW + ^ ffW = 0' (15) 



The first equation implies that cp is an eigenfunction of the Laplace operator on the unit sphere ^ , with eigenvalue 
A^. It is well known that the complete set of eigenfunctions on S^~^ is 



-mD-2 <inD-i < 1710-2, < < mo-s < ■ ■ < mi. 

Here 

3(0) 



Pmi;D-i-ii^) — C*™, ^ (a;) 
is the Gegenbauer polynomial of degree rrii and order (£> — i — l)/2 defined by the generating function [soj : 

1 ^ 

(1 — 2xw + w^) 2 ^.^Q 



and 



The eigenvalue of Ym {0) is 

= 1(1 + D-2), l:^mi. 
For a given / > 0, the degeneracy of the eigenvalue = l{l + D — 2) is equal to 

' - (2/ + i^-2)(/ + i^-3)! 



(D- 2)111 

m2=0m3=0 mu_2=0 mD-i= — 'mD-2 

With A^ = /(Z + D — 2), the equation for g{r) (fT5)) becomes 

AW _^ £+1 ^ _ ^(^ + ^-2)-P-i) ^ 
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whose general solution is given by 



g{r) = r ^ (Ci./,+ ij.^ (wr) + CaiV,^^ (tjr)) , 



(17) 



where J,y{z) and N^{z) are Bessel functions of first and second kind. The third equation in (fTSl) implies that I ^ 
and 



f{r) = 



1 



1 d 



9{r)) = 



1 



1 d 



[r^ CiJi^D^{ujr) + C2Ni^D^{ujr) Y (18) 



l{l + D - 2) dr ^ " l{l + D ~ 2) r^-^ dr 

One can check that g{r) and /(r) given by ([T7| and (fTSl) satisfy the equation ([M)) . As a summary, the set of TM 
modes are given by =0, 

/ D-2 

A, =l{l + D- 2)e— (Ci J,^_z2^K) + C2N,^o^{u;r)) J] ^«";'5i-i(cos^?,) 



3 dr 



(r^ [CiJ,+ £_2(wr) + C2iV,+i2_2K)]) 



i=i 

D-2 



1< ?: < D - 1 



(19) 



It can be written as 2l((/?), where 



CiJ,+i^(^r) +C2A^,+ ^(a.r) J] ^iT;S-.-i(cos^.) e™-- 



is an eigenfunction of the Laplace operator, and 

2t(^), =A^°2-i (r-V) 

a D-2\ dip 



ar 



96*,; 



as defined in [19j, |20|. For fixed ^ > 1, the number of TM modes with nii = I is equal to 5c(/) ([T| 

Next, we turn to explicit solutions of the TE modes. In this case, Ar = and A — e~'"*/(r)A for some transversal 
one form A on the unit sphere. The set of transversal one forms on a unit sphere S^~^ has been constructed in [3ll - [33| 
using more sophisticated language and the results have been employed in (l9l [20| . In fact, we can directly solve the 
system of equations pTt and (IT^ for {Ar,Ag^, l<i<il' — 1}. The form of the equations pT|) and ((T^ allows us 
to solve for Ar, Ag-^, . . . , Aej^_^, A^ one by one using separation of variables. The technique is similar to solving the 
eigenvalue equation for Laplace equations on functions in spherical coordinates. Using the fact that 



de? 



cose, rf^^iTS-.-iCcos^O 



sm tii 



d6. 



+ miinii + D — i — 1) 



mi+i{mi+i + D - i-2) 
siv? 0i 



0, 



one can show that the TE modes can be divided into [D — 2) sets corresponding to {D — 2) polarizations, where for 
1 < j < L) — 3, the i**^ set is given by 

At =Ar = Ae, = . . . = Ae.„i = 

Ae^ ^™.+i(™,:+i +D-^- 2)e— */(r) \{ {su,9,P^^^^^_,icos0,^ _jn,,D^ J Tt 



D-2 



d9i 



n <^i_i(cose,)e' 



i + 1< / < D - 1, 



/(r) = [C^Ji^o^iur) + C^N^^o^itor)) , 

-niD-i < m_D_2 < mo-i, < m_D_2 < • . . < "ii, nii+i > 1, 



(20) 
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and the (D - 2)'*' set is 



D-3 



Pir°".V (cos 0Z,_2) . 

sin 6'd_2 



D-3 



=«e 



i=i 



*/(r) n (sin^,Pi7,^i_i(cos0,))sin0c_2^(p:r.-i(cos0i.-2) 



(21) 



-mo-i < mD-2 < rriD-i, 1 < m£i_2 < • • • < m-i- 



Observe that for 1 < i < D — 3, the i set can be considered as constructed from the transversal one forms on S 
which satisfies Ag. ^ 0. For fixed / > 1, the number of TE modes with m\ = / is given by 



D-i 



where hD.i(l) is the number of modes in the z -set with mi — I, given by 

I mi mi mi+l mo-2 

ho-AD-Y. E- E E - E E 

7712 — 1 — 1 m^_|_i — 1 m^_^2— mo -2—^ rnD-i— — mn-2 
I m2 JTiD-A rtio-s mo -2 

E E- E E El- 

m2 = l m3 = l mD_3 = l 7no-2 = i mo-i= — mD-2 



1, 



hD-D-2(l) 

It is easy to see that 



1 < j < £>- 3, 



/id(0 = Vi(0+E hD-i{m2)^ E (&I^-l(™2) + /^D-l(m2))^ ^ ( ^^"'^ + - 3Z2! ^ + ^i^-i("^2)) 

7712 — 1 7772 — 1 7712 — 1 ^ ^ 

Using this, one can show by induction on Z) and Z that 

l{l + D- 2){2l + D- 2){l + 15-4)! 



/id(0 



{D- 3)1(1 + iy. 



in agreement with the resuh in [3l|. Since we find (D — 2) polarizations for the TE modes, and one polarization for 
TM modes, we have altogether (D — 1) polarizations as expected for electromagnetic field in ZJ-dimensional space. 



Next we turn to the questions of imposing boundary conditions. In the case the electromagnetic field is confined in 
a D-dimensional spherical cavity of radius a, the regularity of the field at the origin r = requires that the coefficient 
C2 in the field modes ((T9| . (|20p and (|2T|) vanishes. The perfectly conducting boundary condition ([7]) is then satisfied 
if and only if for the TE modes. 



J; I D-2 (Uia) = 0, 



and for the TM modes. 



^ / — T / \ 



One can check directly that these coincide with the relative boundary conditions (j!]). On the other hand, the infinitely 
permeable conditions (jS]) are satisfied if and only if for the TE modes. 
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TABLE I: Eigcnmodes of electromagnetic field I: Perfectly conducting boundary condition 



Eigenmodes 



TE modes 
Inside a sphere with radius a 



Lo — ijJ™ , 1,2,... solution of 



(uia) = 



TM modes 



u! — t'-'™, I, j — 1,2, . . . solution of 



/ -D-2 \ 



dr 



= 



Between two spheres with radius a < b 



Eigenmodes uj — uif^ ,1, j — 1,2, . . . solution of 



Degeneracy hoil) = 



J; I D-2 (ljja)N^ ^ D-2 (ijjb) 

- J; I D-2 (Lub)N^ ^ D-2 (iua) = 

1(1 + D- 2){2l + D- 2){l + D - 4)! 
{0-3)1(1 + 1)1 



-{r ^ J,^^(.r)) -(r ^ N^^^(.r)) 



Lu — tiJ;™, I, j — 1,2, . . . solution of 



d I 



d / -D-2 



= 



bo(l) = 



(21 + D -2)(l + D -i)\ 
(D - 2)\l\ 



TABLE II: Eigenmodes of electromagnetic field II: Infinitely permeable boundary condition 



Eigenmodes 



Eigenmodes 



TE modes 
Inside a sphere with radius a 



UJ — I, j = 1,2, . . . solution of 



d / 4-£> 

dr 



(r 2 J^^D^[u:r, 



Between two spheres with radius a < b 



UJ — ujfj^, I, j = 1,2, . . . solution of 



d ( 1^ ^ , 



dr V 

TrV ' 



(u:r) 



d / 

r 2 AT D-2 (tjr; 



Degeneracy 



hD(l) 



1(1 + D- 2)(2l + D- 2)(l + D - 4)! 
(D-3)\(l + l)\ 



TM modes 



UJ — ujij = 1,2,... solution of 



I D-2 (uja) = 



= 



Lj — ciJ;™, I, j — 1,2, . . . solution of 



J, ^ D-2 (uja)N^ ^ D-2 (ujb) 
'J, I D-2 (ujb)Nj ^ D-2 (ujg) = 



bD(l) = 



(21 + D - 2)(l + D - 3)1 
(D - 2)W. 



and for the TM modes, 

J; I D-2 {oja) — 0. 

One can check directly that these coincide with the absolute boundary conditions ([5|). Thus the Casimir energies 
computed in [lOj] can be interpreted as the Casimir energies in the interior of a Z?-dimensional spherical cavity with 
perfectly conducting boundary conditions on the spherical boundary (corresponding to relative boundary conditions) 
or with infinitely permeable boundary conditions on the spherical boundary (corresponding to the absolute boundary 
conditions). 

For electromagnetic fields confined between two concentric spheres with radii a and b, where < a < 6, the 
boundary conditions on the surface r = a and on the surface r = b give rise to a system of equations satisfied by 
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Ci, C2. For perfectly conducting conditions on both shells, we have 



TE modes: 



TM modes: 



dr 



(r^ Ci J; I D-2 (ur) + C^A^; , 0-2 (a;r) ^ 
(r^ [CiJ;+£_2(a;r) + C2Af,+ ii_2K)]) 



= 0, 
= 



The TE modes and TM modes are those uj which give rise to nontrivial solutions of Ci, C2, i-C, they satisfy 
TE modes: J; | 0-2 (uja)Ni ^ 0-2 (ub) — J; | 0-2 {ujh)Ni ^ 0-2 (a;a) = 0, 



TM modes: 



-f 



■D-2 



r—b 



— ir 2 N,,D-2 (ujr 
dr \ 2 ' 



= 0. 



Similarly, for infinitely permeable boundary conditions on both shells, the TE and TM modes are solutions of the 
following equations: 



TE modes: 



TM modes: 



d 
dr 



[r 2 J,^n^{^r)\ -[r ^ N^^n^i^r)] 

r—a 

d ( 4-13 / x\ d ( 4-.P 



= 0, 



The results of this section are summarized in Table|T]and Table HIl Notice that when Z) = 3, h-j,(V) — h-},(V) = 2Z + 1, 
the TE modes for perfectly conducting boundary conditions is the same as the TM modes for infinitely permeable 
boundary conditions, and the TM modes for perfectly conducting boundary conditions is the same as the TE modes 
for infinitely permeable boundary conditions. Therefore when 13 = 3, the Casimir energy for perfectly conducting 
boundary condition is the same as the Casimir energy for infinitely permeable boundary condition. This duality does 
not hold for D > 3. 



III. CASIMIR ENERGY INSIDE A SPHERICAL SHELL 

In this section, we consider the Casimir energy inside a Z3-dimensional spherical cavity of radius a with perfectly 
conducting boundary and with infinitely permeable boundary. As mentioned in the previous section, the Casimir en- 
ergy has been computed in [l9l | under the context of absolute and mixed boundary conditions using zeta regularization 
method. Here we compute the Casimir energy using cut-off regularization. 

Using the cut-off regularization method, the Casimir energy inside a spherical shell of radius a is given by 

where A is a cut-off parameter. Introducing the zeta function 

C\s) = Qf^^s) + C^Us) = J2 W E + E ^^(0 E (^™)"'' ' (22) 

1=1 j=l 1=1 3 = 1 

it is standard to show that up to the A*' term, 

pint . N ^^ r(£> + l-ii) int,»-n-l V'(l)-l0gA i„t , It^p Aint. n 

-E^Cas(a) = 2^ Yl^D-^) ^^+1 ^ 2^=-5^ ^ 



11 



TABLE TIL Heat kernel coefficients for tlie interior 



of a splierical shell of radius a with perfectly conducting boundary 
condition 





D = 3 


D = A 


D = 5 


D = 6 


D = 7 


D = 8 


D = 9 


int 

Co 




3a* 


a^ 


5a'' 


a' 


7a« 


a» 


30F 


32 


150F 


384 


14O0F 


6144 


189O0F 







^a^ 


a* 


3v7Fa^ 


a" 


5^o7 


aS 


16 


12 


128 


60 


1536 


560 


„int 
C2 


4a 


3a2 


2a^ 


5a* 





7a'' 


a' 


3^^F 


8~ 


9^ 


192 


2304 


315^ 




5 

8 


21l7¥a 
512 


a^ 
2 


5850Fa^ 
4096 


5a* 
48 


lO150Fa* 
49152 


190" 
1680 


„int 
C4 


16 


7 


796a 


3a^ 


6080=^ 


197a* 


346a5 


315v'7i"a 


20 


945 v/¥ 


8~ 


15750F 


1920 


47250F 


int 
C5 




16310r 
65536a 


541 
2880 


75361-^/7Fa 
524288 


163a2 
768 


6994813^™^ 
94371840 


887a* 
13440 


int 






6632 


1109 


17468a 


a' 


1915624a^ 


C6 






45O450Fa 


7560 


614250F 


~ 8 


141891750F 


C7 








1052991^7? 
16777216a 


143263 
967680 


23185O1770ra 
3019898880 


491a^ 
5120 


int 
C8 










16438144 
727476757¥a 


33521 
226800 


-0.12483a 


„int 

Cg 












2580976217486942940567943799^^ 
32674585544991625157478973440a 


0.14319 



0.15150 



a 



where 

C*=Res,^i2^ (r(s)r*(s)), (23) 

and FPg^_ 1 C'"'(s) is the finite part of the zeta function ^"'(•s) at s = — -. In the zeta regularization scheme, the 
terms with negative powers in A are omitted. The zeta regularized Casimir energy is defined as 



cifin (I . ,\ 1. 



pint, zeta/ \ _ M ,int ( ^ 

1 



2 ^ V 2 

Here /i is a normalization constant. The 1/e + log/x^ term is the ambiguity arises when C'"'('5) has pole at s 



Systematic ™.h„d ha, b»n developed i„ MUMa tor eotnputing the zeta f„„eti„» i„ spheriea, ge„„etriel 

Since the computation is quite involved, we leave it to Appendix [X] 

Using the fact that the Hurwitz zeta function C,h{s\~x) (lASp has only one pole at s = 1 with residue 1, we can 
readily find the heat kernel coefficients ([23l) from the expression for the zeta function C'^H'S) (|A9p . The results for 

3<D<9,0<i<L' + lare hsted in Table |IlT] and Table HV] The finite part of the zeta function C'"'(s) at s = -- 

can also be calculated directly from (|A9[) . The result is listed in Table |Vl For 3 < £) < 6, we find that our result is in 
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TABLE IV: Heat kernel coefficients for the interior 



of a spherical shell of radius a with infinitely permeable boundary 
condition 





D = 3 


_D = 4 


D = 5 


D = 6 


D = 7 


D = 8 


D = 9 


int 

Co 




3a* 


a^ 


5a'' 


a' 


7a« 


a« 




32 


15V? 


384 


1407? 


6144 


189O0F 







V^a^ 


a* 


30Fa5 




5-y7ra'^ 


a» 


16 


12 


128 


60 


1536 


560 


„int 
C2 


4a 


3a2 


2a^ 


5a'' 





7a« 


a^ 




8 


9x/;f 


^192 


2304 


3157? 




5 
8 


121^^ 
512 " 


4 


235 y^a^ 
4096 


7a* 
240 


49152 


a« 
1680 


Jnl 

C4 


16 


49 


68a 




16a^ 


9a* 


82a'5 


SlS-y/Tra 


60 


189v/? 


8 


175^/? 


640 


33O750F 






27130r 
65536a 


2429 
2880 


44071 0ra 
524288 


23a^ 
256 


2O365870ra^ 
94371840 


51a* 
4480 


ce 






11048 


6199 


58012a 


2 

a 


90589oa 






45045 0Fa 


7560 


2252250F 


^12 


141891757? 










8713390F 
8388608a 


785567 
967680 


287752910Fa 
377487360 


75a2 
1024 












2955168 


185449 


-0.12063a 










8083075V¥a 


226800 



4"* 



int 
ClO 



4030878578159017023021286121v/7r 
' 32674585544991625157478973440a 



0.82406 
0.22939 



good agreement with the result in ^J^. Notice that ^"^(s) = a^*C™*(s), where C'"*(s) is independent of a. Therefore, 



FP,=_i/2r*(s) = FP,^_i/2r'(s) +21ogaRes,=„i/2r*(s) = FP,^_i/2r*( 



int 



logo. 



The presence of the log a term indicates the ambiguity in defining zeta regularized Casimir energy. 



IV. CASIMIR ENERGY IN THE ANNULAR REGION BETWEEN SPHERICAL SHELLS 



In this section, we consider the Casimir energy inside the annular region between two spherical shells in D- 
dimensional space with radius a and 6, where a < b. The Casimir energy outside a Z?-dimensional spherical cavity 
with radius a is obtained by taking the limit 6 — >■ 00. 

As in the previous section, the Casimir energy inside the annular region between two spherical shells with radius a 



13 



TABLE V: The value of FPs^_i/2C'"'(s) for the interior of a spherical shell of radius a with perfectly conducting and 

infinitely permeable boundary conditions 



perfectly conducting 



infinitely permeable 



D = 3 



_D = 4 



D = 5 



D = 6 



D = 7 



D = 8 



D = 9 



0.16785 _^ 16 log g 

a 3157ra 

0.04881 1631 log a 
a 65536a 

0.01880 6632 log g 
a 450 457ra 

0.02022 1052991 log g 
a 116777216g 

0.01603 16438144 log a 
a ~^ 727 476 757ra 

0.00591 _ 2580976217486942940567943799 log a 
a 32674585544991625157478973440g 

0.00358 0.08547 log a 
+ — 



0.16785 _^ 16 log g 
a 3157rg 



0.34619 ^ 2713 log g 
g 65536a 



0.46773 11048 log a 
a 450457rg 



0.52494 871339 log a 
a ^ 8388608g 



0.56057 2955168 log a 
a 80830757rg 

0.59822 4030878578159017023021286121 logg 
a ^ 32674585544991625157478973440a 



0.63741 0.12942 log a 
+ — 



and b, where a < &, is given by 

OO CXD OO OO 



1=1 

D-l 

E 



1=1 ]=i 



rp + l-z)^,„„^,_^_i V(i)-iogA 



' anil \ 

-Ca a 



J ^^to ' ' anil I 1^1) /^ann / „\ 

^^7^c^+i + 2FP.=-4C is). 



As before, 

c^"" = Res^^z2_, (r(s)C""(s)) . 
Using the same approach as in Appendix [Xl we find that for perfectly conducting boundary condition, 

OO OO 



1=1 



for infinitely permeable boundary condition, 



(24) 



where 



1=1 



1=1 



(25) 



^log^[cJ^{az) + azJl{az)] [cN^{bz) + hzNl{bz)] - [cj^{hz) + bzJl{hz)] [cN^{az) + azNl{az)]^ 



dz. 
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Using the fact that Hi^\z) = J^{z) + iN^{z), H'if' {z) = J^{z) - iN^{z), J^iiz) = e"^ ^{z), Jy{~iz) = 

Ii,(z), Hl^\iz) — — K^(z), H^-^h—iz) = e^K^(z), where /^(z) and K„(z) are the modified Bessel 

functions of the first kind and second kind (see [30|), one finds that 



.3,nn,i., ^ ,. sin(7rs) 

sin(7rs) 



m— >0 TT 

d 



(z^ - m^) "-^ log^Ii,{bz)K^{az) - I^{az)K^{bz)^dz, 



X 



log^[cI^{az) + aziliaz)] [cK^{bz) + bzKl{bz)] - [cl,y{bz) + bzll{bz)] [cK^{az) + azKl{az)]^dz, 



for - < Re s < 1 . Notice that these can be written as 
2 



>ann , i/ / \ >int .v / \ , >cxt . i/ / \ , /-inter , i/ / \ 

Cd (■s) = Ci3 (s)+Cd ■ {s) + Cd (s), 

j-a.nn.i' / \ Ant.i' / \ , >cxt.i/ / \ , /-inter, i/ / \ 

Cr,c is) = Ci?,c (S) + Ci?,c (■«) + Cr,c (■«), 



where 



sin(7rs 



m— >0 TT 



Cr'''(.)= hm^ 



sin(7rs 



/-inter, iv" / \ i • 



m— >0 TT 

sin(7rs 



?7i— >-0 TT 

sin(7rs 



m— J-O TT 

sin(7rs 



m— >0 TT 

sin(7rs 



oo , 
2 ™2^-s " 



{z'-m'y 



dz 



logjz-''/, (6z)}dz, 



°° d r ^ 

2 ™2^-s log|z''X^(az)}dz, 



dz 



oo 



iz'-m')-'^iog^l- 



I^{az)K^{bz) 



dz, 



oo , 

dz 



>dz. 



m— >0 TT 



/j,(6z)ii',,(az) 
{z-''[c^(&2) + 6z/;(6z)]} 

(z^ - m^)"''-^ logl^" [-cK^{az) ~ azKl{az)]^dz, 

[cl^{az) + azll{az)] [cK^{bz) + bzKl{bz)] 



oo , 
2 ™2n-s " 



{z^ ^m^y^lofii 1 
az 



[c/y(6z) + bzll{bz)] [cKjj{az) + aziir^(az)] 



dz. 



Correspondingly, the zeta function ^'"^(s) (|24l) can be decomposed into three terms 

which gives rise to the decomposition of the Casimir energy: 

i?S:^(a, b) = E^IM + + E^^\a, b). 

As seen from Appendix \K[ E™}^^{b) is the Casimir energy inside a spherical shell of radius b. When & ^ oo, this can 
be considered as the Casimir energy of the free space. On the other hand, we show in Appendix |A] that in the limit 
b — oo, 

hm CT'^'^is) = hm Ck-^f'-'^Cs) = 0. 



Therefore, E'^^^{a) is identified with the Casimir energy outside a spherical shell of radius a (which has been renor- 
malized by subtracting away the Casimir energy of the free space). We call i5Q*g'^(a, 6) the interacting term of the 
Casimir energy between the spherical shells. 

We consider the Casimir energy outside a spherical shell £^q^* (a) and the interacting term EQll''{a, b) separately in 
the following sections. 



V. CASIMIR ENERGY OUTSIDE A SPHERICAL SHELL 



The computation of the Casimir energy outside a spherical cavity of radius a follows the same way as the computation 
for the Casimir energy inside the spherical cavity. We leave it to Appendix [Xj It can be shown explicitly that for 
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TABLE VI: Heat kernel coefficients for the exterior of a spherical shell of radius a with perfectly conducting boundary 

condition 





D = 3 


D = 4 


_D = 5 


D = 6 


D = 7 


D = 8 


D = 9 


int 


4a 


211^7ra 


796a 


7536iy7ra 


17468a 


231850177 V7ra 


0.12483a 




512 


945 


524288 


614250F 


3019898880 


cxt 
Cd 


3 

^8 


13 

~20 


2339 
2880 


6451 
7560 


824417 
967680 


193279 
226800 


-0.85681 


int 


16 


163iy¥ 


6632 


1O529910F 


16438144 


25809762174869429405679437997^ 


0.15150 


315v'7ra 


65536a 


45045 0Fa 


16777216a 


727476750Fa 


32674585544991625157478973440a 


a 



TABLE VII: Heat kernel coefficients for the exterior of a spherical shell of radius a with infinitely permeable boundary 

condition 



D=3 D=A D=5 D=Q D=7 D=S D=9 



4a 


121v^a 


68a 


44071 V^Fa 


58012a 


28775291^™ 


0.12063a 


30F 


512 


1897¥ 


524288 


2252250F 


377487360 


3 


11 


451 


1361 


182113 


41351 


-0.17594 














8 


60 


2880 


7560 


967680 


226800 


16 


2713^^ 


11048 


8713390F 


2955168 


40308785781590170230212861217¥ 


0.22939 


315\/7ra 


65536a 


45045^/¥a 


8388608a 


8083075V^a 


32674585544991625157478973440a 


a 



< i < D — 2 and for i — D — 1,D + 1 and D odd, the heat kernel coefficients for the exterior of the sphere c^''* is 
related to the heat kernel coefficients for the interior of the sphere cj"* by 

c^^^t = (_l)*+icj"'. (26) 

For i = Dori = D — 1,D + 1 and D even, this is not necessary true. The values of c™' for 3 < £> < 9 and 
j = D - 1, £>,£> + 1, is listed in Table EH and Table |VlIl Compare to Table Im] and Table IVl we find that ^ does 
not hold ioi i — D, but it still holds i — D — 1, D + 1 and D even when 3 < Z? < 9. It seems that the latter will 
be true for all D > 3 but we don't know a proof for it when D is even. 

The value of the finite part of the zeta function C'^'^'(s) at s = — - can be computed from (|A13p and the results for 
3 < L» < 9 are listed in Table EnH 



VI. TOTAL CASIMIR ENERGY OF A SPHERICAL SHELL 



In this section, we consider the Casimir energy of a single spherical shell system. The Casimir energy for a perfectly 
conducting spherical shell or an infinitely permeable spherical shell is the sum of the Casimir energy inside the spherical 
shell and the Casimir energy outside the spherical shell. In the cut-off scheme, we find that 

psinglc/ N _ T (D + 1- i) single .^i-g-l '0(1) - fog A single , 1 ^single/ \ 

^Cas W - 2^ T[^-') * + 2^^s=~^^ [S), 

up to the term constant in A. Here 

^single ^ ^int ^ ^ext ^ ^single ^ ^i„t ^ ^ext 

(1^5)) shows that if i is even and i ^ D, c^'"s'° — 0, but if i is odd, c™^'° is in general nonzero. However when 
D = 3, c^^^^° = since both c'"' and vanish. Therefore, when D = 3, Cq"^^'', c^"®'°, €2"®'°, €4"®'° all vanish and 
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TABLE VIII: The value of FPs=_i/2C'''''(s) for the exterior of a spherical shell of radius a with perfectly conducting and 

infinitely permeable boundary conditions 



perfectly conducting 



infinitely permeable 



0.07549 _ 16 logo 
a 3157ro 

0.18230 _ 16311oga 
a 65536 

0.29414 _ 6632 logo 
a 45045770 

0.38128 _ 1052991 logo 
a 1167772160 

0.44667 _ 16438144 logo 
a 727 476 757ra 

0.50087 _ 2580976217486942940567943799 log a 
a 32674585544991625157478973440a 

0.54880 0.08547 log a 



0.07549 _ 16 logo 
a 3157ro 

0.02446 _^ 2713 logo 
a 65536a 

0.01042 _ 11048 logo 
a 450457ro 

0.00447 871339 logo 
o 8388608O 

0.00566 _ 2955168 logo 

a 80830757ra 

0.01784 4030878578159017023021286121 log o 
o 32674585544991625157478973440a 

0.02920 0.12942 logo 



TABLE IX: The value of — FPs=_i/2C'*'°^'''(s) for a spherical shell of radius a with perfectly conducting and infinitely 

permeable boundary conditions 



perfectly conducting 



infinitely permeable 



0.04618 



0.11555 _ 1631 logo 
a 65536 

0.13767 



o 



0.20075 _ 1052991 log a 
a 116777216a 

0.21532 



a 



0.04618 

a 

0.18533 _^ 2713 log a 
o 65536a 

0.22866 



a 

0.26470 871339 log a 
a 8388608a 

0.28312 



0.25339 _ 2580976217486942940567943799 log o 
o 32674585544991625157478973440a 

0.27619 



a 

0.29019 4030878578159017023021286121 log o 
a 32674585544991625157478973440a 

0.33331 



one can take A — )• 0+ in the cut off regularized Casimir energy to obtain an unambiguous Casimir energy given by 

-FPg=_i/2C^"^^^'^(*) whose numerical value is well-known to be correct to five decimal places. When D > A, 

divergence is always present in the cut-ofF regularized Casimir energy. The leading divergence is of order A~^. In 
order to obtain a physically meaningful vacuum energy, one need to remove the divergence by some renormalization 
procedures. 
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In the zeta regularized scheme, the Casimir energy is defined as 



single 



1 . o\ I 



The value of -FPs=_i/2C'^'"^'''(s) for 3 < Z? < 9 is tabulated in Table ITXl As mentioned in the previous section, when 
D is odd, it can be shown that c^^^ — — c'^^+i- Therefore, c^^f^i ~ when D is odd. In this case, the zeta regularized 



The Casimir force acting on a single spherical shell is given by 



Casimir energy is finite and given by -FPs^_i/2C"^''^('S)- 



Using cut-ofF regularization, we find that there are divergences in the Casimir force for D > 4. Using zeta regularization 
scheme, there is no divergence when D is odd, and the Casimir force is given by 

T-i single, zeta/ \ 
^single, zeta j-^-j _ -^Cas 

When Z3 = 3, we have the well-known result that the Casimir force is repulsive. However, for 4 < D < 9 an odd 
dimension, we notice from Table IIXI that the Casimir force is attractive for perfectly conducting boundary condition, 
and repulsive for infinite permeable boundary condition. 



VII. CASIMIR EFFECT OF DOUBLE SPHERICAL SHELL CONFIGURATION 



In this section, we consider two concentric spherical shells of radius a and 6 respectively, where a < b. The total 
Casimir energy of such a system is the sum of the Casimir energy inside the sphere of radius a, the Casimir energy 
inside the annular region bounded by the two spheres, and the Casimir energy outside the sphere of radius b. Namely, 

i?c°as"°(«, b) = E'SM + E^^J^ia, b) + E^tib). 
As discussed in Section ITVl the Casimir energy in the annular region i?™g(a, b) can be decomposed into 

E^l^ia, b) = E'SUb) + i?cas(«) + E'SlT'i^, b). 
Therefore, we can rewrite the total Casimir energy of a double spherical shell system as 

£^g°af' ^(«, b) = E^:f\a) + El^tria, b) + E^^^ib), 

which is the sum of the Casimir energy of a single spherical shell of radius a, the Casimir energy of a single spherical 
shell of radius b, and the interacting term E^(^^°'-'{a, b). The Casimir energy of a single spherical shell has been discussed 
in the previous section. We find that using cut-ofF regularization, it does not have divergence only when Z? = 3. Using 
zeta regularization, the Casimir energy does not have divergence for all odd D. When D = 3, both the cut-off 
regularization and zeta regularization yield the same result. 
Using cut-off regularization, we have as before, 

up to the term constant in A. Here 

^double = C''"^'"(s; a) + C"Slc(s. ^) ^ (into'(s). 

In Appendix [XI we show that C'"*'^''(s) is analytic for all s. Therefore, when 0<i<D— loi-i = D + l, the 
contribution to 

^double ^ ReS^^l.-. (r(s)C'^°""'=(s)) 
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comes only from C'''"sie||g. gj^^ ^smgiej-g. f^y Therefore, the A — )■ 0+ divergence term of the Casimir energy £'cas'^'°(^) 
is a sum of the divergence term for a single spherical shell of radius a and the divergence term for a single spherical 
shell of radius b. As a result, the divergence vanishes only if = 3. 

Since the single shell contributions to the Casimir energy have been discussed in the previous section, we focus our 
attention on the interacting term. From the result of (jA17p in Appendix [Xj we find that for perfectly conducting 
boundary condition, the interacting term of the Casimir energy is 



1 ... 



dz 



1 = 1 -"^ 



logU 



(az) + azi' {az) ^K,^i) (bz) + bzK' (bz) 



For infinitely permeable boundary condition. 



E'S::^{a,b)^-J2hn{l) I log<jl 



1=1 



(bz) + bzl' {bz) i^K.^i) {az) + azK' {az) 



{az) + azr^,^ {az) ^K^n) {bz) + bzK^^^^ {bz) 



(27) 



dz. 



dz 



—h{i){bz) + 6zr(,j(6z)J ^-^K^^i){az) + azKl^^^{az) 

^ 1 V^^ m ^^ /i lHi)i'^^)K.(i){bz) \ ^ 

+ ^1^^d{1) log<^l-- — — i—r>dz. 

27r^ Jo { I^(i){bz)K^^i){az) j 

The contribution of the interacting term to the Casimir force acting on the shell with radius a is given by 

d 



(28) 



Cas 



'■(a, 6) 



da 



In case of perfectly conducting boundary condition, 



FclTi'^^b)=—J2hi{D)J^ A''^'+^{z;a,b)dz+—J2biiD)j^ A^^'lZf {z; a,b)dz. 



In case of infinitely permeable boundary condition, 

F^'^"^^{a^b)^—J2hi{D)J^ A';[t^{z;a,b)dz+—J2bi{D)J^ 



a,l+- 



(z; a, b)dz. 



1=1 



Here 



A 



a,u{l) 



A 



R,c 



{z;a,b) 



{z;a,b) 



Ku(i){bz) / (i^ h(i){az)K^(i){bz) 
Iu(i){bz)K^(i){az)'^ / \ I^(i){bz)K^(i){az) 

_ -[c-R-,(i)(bz)+fa<(,)(bz)](a"^" + (t'"-c")) / 



^ _ Y-I.^(l){o.^)+o-^ll{l){o.z)\[':K^(i){bz)+bzK'^^i-^(bz)\ 

[c/i,(;) (fcz)+fczr (fcz)J ^cK^,^i) (az)+azKl^^^ (az)J J 



u{l){'^>z)+bzIl^^^{bz)\^ ^cK„^lS|(az)+azKl^^^^(az)^ j 

Notice that the function 

ch(i){z) + zll^i^{z) = {c + v{l))I^^i){z) + z/^(,)+i(z) 

is > when z > 0. For c — ^^-^ or and v{l) = I + ^{l) — c is always positive. Therefore the function 

cK^(^i){z) + zKl^i^{z) = (c - v{l))K^i^i){z) - zK^^i)_i{z) 
is < when z > 0. On the other hand, since 

d K^(i){z) 1 



dz Iv(i){z) zl^i^i){z)'' 



<0, 



d 

dz 



cK^(j){z) + zK',iAz) 



z^ + {u{l)^ - c2) 



c4(,(z) + z/:(,(z) ; ,(e/,(,(,) + ,/^^^^(,))^ 



< 0, 
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Iiy(i)iz) and — (^cKi,(^iSi{z) + zKl(^i-^{z)^ ^ (^cl^(^i){z) + zllf^i^^lz)^ are nonnegative decreasing functions of z 
when z > 0. Therefore, when z >0, v{l) = I + ^^-^ and c = or we find that 

^ ^ K^(i){bz) ^ K^(i)(az) ^ ^ cK^ji) jbz) + bzKl^^^ jbz) ^ cK^(i) jaz) + azKl^i^ jaz) 

Iu{i){bz) ~ I^(^i){az)' ~ cIi,Q^{bz) + bzll^i.^{bz) ~ cl^i^i){az) + azll^^^^{az) 

These imply that A^'^'''^(z; a, 6) and A^'^'''' (z; a, 6) are nonnegative functions of z when z > 0. As a result, we find 
that for either perfectly conducting or infinitely permeable boundary conditions, the interacting term of the Casimir 
energy gives rise to a force acting on the inner shell that tend to push the inner shell towards the outer shell. 

Next we consider the interacting term of the Casimir force acting on the outer shell with radius b. This is given by 

For perfectly conducting boundary condition and infinitely permeable boundary condition, we have respectively 

PcT'(^^b)=-—S2hi{D) ^'^+^^^.a,b)dz-—Y,bi{D) A'';;ff{z;a,b)dz. (30) 
1=1 •'^ 1=1 •^o ' ^ 

and 



Here 



Fc'r\<^,b)=-—l^hiiD) A"^ll^{z;a,b)dz-—l^kiD) I A^^^ (z; a, (31 



1 ^ poo 



.b,u{i)f ,^ Iu{i){az) I ( Iy{i)(az)K^(i){})z) 

^ ^ X,(,)(az)4(;)(6z)2/ 1^ I,^l){bz)K,^l){az) , 

It follows from (1^^ that A^"^''-* (z; a, 6) and A^'"^*-'-* (z; a, 6) are always nonnegative when z > 0, = / + and 

c = -^^^ or ^-j^. Therefore, the minus signs in front of the integrals in (l30l) and ((3T|) imply that the interacting term 
of the Casimir energy always gives rise to a force acting on the outer shell that tends to push the outer shell towards 
the inner shell. Combining with the result about the interacting term of the Casimir force on the inner shell, we 
find that the interacting term of the Casimir energy gives rise to a Casimir force that attracts the two shells to each 
other. Taking into account the single shell contribution, the Casimir force can only be unambiguously regularized in 
the context of cut-off regularization when the space dimension is Z? = 3. In this case, we see that acting on the inner 
shell with radius a, the total Casimir force is 

j-ia [ i\ 77isine;le/ \ . T-^a, inter/ t\ 0.04618 ^-^a. inter/ 7\ 

^Cas(a, b) = J (a) + F^^^ (a, b) = — — + F^^, (a, &), 
Therefore it is always positive (pointing outward). However, for the outer shell of radius 5, the total Casimir force is 

771b / T\ 7-isinglc/i\ I 771b. inter/ i\ 0.04618 T-^bintcr/ i\ 

Fce.s {a,b)^ F^J (6) + F^^^ {a,b) = — + F^'^^ (a, 6) . 

The first term is positive but the second term is negative. It is easy to show that when 6 — > A^"*"'' (z; a, 6) and 
A^'^J'-' (z; a, &) approach positive infinity, and when b — > oo, they approach zero exponentially fast. Therefore, for b 
close to a, the Casimir force acting on the outer shell FQ^^{a,b) is negative (pointing inward). However, for fixed a, 
when b is large enough, the Casimir force will become positive. There is an equilibrium point where the force is zero 
but it is unstable. Numerical simulation shows that the equilibrium point appears at b/a = 5.55857 correct to five 
decimal places. If the space dimension D is odd, one can also discuss the sign of the total Casimir force acting on the 
two shells in the similar way under the context of zeta regularization. 
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VIII. ASYMPTOTIC BEHAVIOR OF THE CASIMIR FORCE FOR SMALL SEPARATION OF SHELLS 

In this section, we consider the asymptotic behavior of the Casimir force when the separation between the plates 
d = b — a is small compared to the radii of the spheres. Since the single shell contribution to the Casimir force does 
not depend on the plate separation, we only consider the interacting term. We first consider the force acting on the 
inner shell. By a change of variable, we find that in case of perfectly conducting boundary conditions, the interacting 
term is given by 



a ' ' a ' 



(32) 

Using Debye uniform asymptotic expansions of the Bessel functions (|A4[I and (|A11[) . we find that as v ^ oo, 



n—1 \ i—1 



(33) 



Pni{z; a,b) and qni{z;a,b) are polynomials of t{z) = l/Vl + z'^ and t{bz/a) that vanish when b a. Using the 
convention Pnoiz; a, 6) = qno{z', a, 6) = 1, ([33]) implies that 



1 D-2 / \ 

™ 1 = 1 n=l i=0 ^ 



—. dz 

* J ti-isx> \ \ ^ / / ^ 

= — — / r(.s)2-'^ EE E + « - J - 2; f ) / ( ( - j - r7(^) ) Pm{z; a, b)dzds 

(34) 



for large enough /i. Notice that 





fhz 












fbz 


V { 


K a 



ri{z) — log — as 2; ^ 



a 

rjiz) ^ ( 1 I 2 as z — oo. 

\a J 

Therefore the integral in z in the last line of ([34]) is convergent as long as s > 2. Let 

_ 6 ^_d 

The leading term of p4p when 7 — > is given by the term with i — Q,j = D — 2 which gives rise to a pole of 
Ch {s + i — j — 2; y) Sit s ~ D + 1. We find that this leading term is 

1 T{D +_1) g ^-n-i^^^^_^ r f fbz\ ^^^A -^-^ ^3^^ 



7ra2 2^T- Ja \ \ a 

n—1 ^ 



By definition xd-j — — 3)!. On the other hand, to the first order in 7, 

V (^^) - viz) = - ^) '?'(^) + o (7') = iVTT^ + o (7') 
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Therefore, the leading term of (1351) when 7 — is given by 



ira^jD+i 2D+1 '{D- 3)1 Jo (l + z2)* V^a^^^+i 3^+1 {0-3)1 T (f ) ' 

In the similar way, one can show that the leading contribution of the second term in (j32p when 7 — is given by 

1 r{D + 1) Cb.{d + 1) r (^) 
^^2^15+1 2^+1 iD-2)i r(f) ■ 

Summing up, we find that as the separation between the shells d goes to zero, the leading term of the Casimir force is 



,a.inter, 1 T {D + 1) UD + 1) T {^) , 1 T {D + 1) + 1) T ( ^'^^ 



pa, inter/ t\ ^ ^ \^ ' ^/ ^ ^ ^ a. ^ ' ^^^^ ' ^/ \ 

Cas ^^a2^D+i 2^+1 (I? -3)! r(f) V^a27^+i 2^+i (7^-2)! T 

A{S"-')DT{^)CRiD + l) 



2 
> 2 ; 



(36) 



= (i^-l 
where 



-D + l 



D- 



2 , 



is the surface area of the sphere of radius a. For infinitely permeable boundary condition, the same derivation shows 
that the leading term of the Casimir force when 7 is small is also given by ([55)) . If we consider the Casimir force 
acting on the outer spherical shell, then the leading term is give n by the negative of (1361) . which signifies a force in 
the opposite direction. It is well-known that (see e.g. Q and |39|) 

OT(^)Cfl(i^ + l) 1 

is the magnitude of the attractive Casimir force per unit area acting between two parallel plates due to massless scalar 
field with Dirichlet boundary condition. Multiplying by the factor {D — 1) gives the magnitude of the Casimir force 
density due to electromagnetic field, where the factor [D — 1) accounts for the different polarizations of photon in D 
dimensional space. Therefore the analysis in this section shows that in the limit 7 = {b~ a)/a is small, one finds that 
the leading term of the Casimir force per unit area between two concentric spherical shells is the Casimir force per 
unit area acting on a pair of parallel plates. This is a verification of the proximity force approximation [40l l4l|. 

The lower order corrections to the proximity force approximation can be computed from p4p by expanding V {^) ~ 
r]{z) and Pni{z] a, b) in terms of 7. The computation is quite involved and we leave it to the interested readers. 



IX. CONCLUSION 



This work considers electromagnetic field in Z3-dimensional spherically symmetric cavities subject to perfectly 
conducting or infinitely permeable boundary conditions. We derive explicitly the eigenmodes of electromagnetic field 
in Z?-dimensional spherical coordinates. The field modes are divided into TE modes and TM modes, and it is shown 
that TE modes can be decomposed into {D — 2) polarizations, whereas TM modes only have one polarization. In 
total, there are (D — 1) polarizations as one should expect. Although perfectly conducting and infinitely permeable 
boundary conditions are defined in terms of the two form F^^dx^dx'^ which represents the strength of the field, 
it is shown that in D-dimensional spherically symmetric cavities, the perfectly conducting and infinitely permeable 
boundary conditions are equivalent respectively to the relative and absolute boundary conditions on the potential one 
form A = A^.dxf' defined by Gilkey [201. 

The Casimir energies in the inside and outside of a spherical shell are computed using both cut-off regularization 
and zeta regularization. The sum of the inside and outside contributions is the Casimir energy of a single spherical 
shell system. It is observed that when the space dimension is D = 3, the Casimir energy computed using either cut-off 
regularization or zeta regularization are divergence free and hence agree. When > 4, divergences always present 
in cut-off regularization. However, using zeta regularization, divergences are only present when D is even, as is the 
case for massless scalar field. The heat kernel coefficients that appear in the divergences in the cut-off regularized 
Casimir energy are tabulated for 3 < I? < 9. Casimir energies of the interior and exterior, and the Casimir energy 
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of a single shell are also tabulated for 3 < D < 9. The results disagree with those tabulated in [l|, 0, for the 
case of perfectly conducting boundary conditions since the TE modes are equated to modes of massless scalar field 
with Dirichlet boundary conditions in P, [3 H^, which does not give rise to the correct number of polarizations. 
However, when 3 < D < 7, the values of the zeta regularized Casimir energies inside the spherical shell agrees with 
those computed in [19]. 

For the system of two spherical shells with radii a and 6, the total Casimir energy can be written as a sum of the 
single shell contribution from the inner shell, the single shell contribution from the outer shell, and an interacting 
term. It is shown that the interacting term always gives rise to a Casimir force that attracts the two shells to each 
other. When the separation between the shells d — h ~ a \s much smaller than the radii of the shells, the Casimir 
force behaves like {d/a)~^~^ . The leading term of the Casimir force per unit area is equal to the Casimir force per 
unit area acting between two parallel plates, as what one would expect from proximity force approximation. 



Appendix A: Zeta Functions 

In this section, we conipute the zeta functions that appear in the expressions of the Casimir energies. We follow 
the same approach as in [l^, [l^, 133436} . 



1. Zeta function for Casimir energy inside a spherical shell 



For Casimir energy inside a spherical shell, the corresponding zeta function is 

r\s) = c^FUs) + CiFLis) = E hoii) E «)"'^ + E boil) E (-™)"'' ' (Ai) 

1=1 J=l 1=1 j=l 

where and i^j™ are the TE and TM modes (see Table U and Table HI)) . Using residue theorem, we find that for 
perfectly conducting boundary condition. 



1=1 1=1 ' ^ 

and for infinitely permeable boundary condition. 



c™(s) = E'^^(oC)!if^(s), c™(s) = E^^(^)cr'^'^(s)- (A3) 



J.int,^^ 
Od(1)C 

1=1 ' 1=1 



Here 

>dz 



where 7 is the contour consists of the line a: + ie, .t = cx) — > 0, followed by the semicircle ee , ^ ^ 2' ^ "2"' ^'^^^^^^'^ 
by the line x — ie,x — ^ 00. The factor z^"^ is introduced to cancel the zeros of Ji/(z) at z = 0. The constant 
m > is introduced to simplify the manipulation of the integral. At the end, m is set to zero. Using the fact that 
Ji,{iz) — e^Ii,{z), J^{—iz) = I„{z), where I,j{z) is the modified Bessel functions of the first kind (see [s^]), one 
finds that 



G*"(^)= lim r\z'-m')-^^log{z-''lAaz)} 



dz 



Q'^'^is) = hm / ^\z' - m2)-^-f loglz"^ [c/,(az) + aziUaz)]] 



dz 



smlTTs) / I 
lim ' ' ' 



m->0 TT /am, \ \ a 



1 

— m 



^logjz " [clt,{vz) + i^zll{i'z)]^dz, 
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for — < Res < 1. To obtain the analytic continuation of '"(s) and CRc'^i^)^ one needs the Debye uniform 
asymptotic expansions of the Bessel function Ii,{vz) and Il{vz) (37l Issj: For fixed z, as — >■ oo, 

1 e'-'(^) ^ uu{t{z)) 



/2^ (1 + ^2)1/4 II 
k—0 



12-EV 



k=0 



where 



T]{z) = + Z^ + log 



1 + 



t{z) 



and Ufc (t) and Ufc (t) are defined recursively by 



Uoit) = 1, life (t) - + - 



(1 - 5T^)ufe_l(T)dT, 



clv{vz) + vzllivz) 



/^e''''W(l + z2)i/4 



27r 



ct(2;)ufe_i(i(z)) + Vk{t{z)) 



Let 



log 



>^ Mfe(t) 



fc=l 



E 



loe 



ctUk-l{t) + i;fe(t) 



fc=l 



E 



M„,c(i) 



It can be shown that Dn{t) and Mn{t) are polynomials of i of the form: 



(A4) 



Vo{t) = 1, Vk{t) = Uk{t) - t^l - t^)Uk_^{t) - ^ ^ ^ Mfc-lW- 

It can be shown that Uk{t) and i'fc(i) are polynomials in t of the form ajt^^'^K For c/i,(i/z) + i^z/^(i/z), (jA4p implies 
that 



(A5) 



M„(t) = ^m„j(c)r+2J'. 



From (|A4|) and (IA5|) . we have 



-log{z-'^/.(^.z)} 



z dz 



■ log (e-(^)) - ^ log(l + zY^ +^^J2 + i (^ogU.z) log 



4=1 



/2^(l + z2)l/4 



N 

E 



A(i(^)) 



1 z 



N 



2 1 + ^2 



E ,; + - ( ^ogUi^z) - log ■ 



/2^(l + z2)l/4 



E 



A(t(^)) 



(A6) 



Similarly, 



1 z 



dz 



\og[z~'' [ch{yz) + yzllivz)]] (^/l + z^- l) + 3 j 



ztjzfMlMz)) 

y2 



I/' 



+ ^ (log [clu{yz) + vzll{vz)] - log ^ >- Yl y 



(A7) 
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Here N is any positive integer. Using (jA6p and (jA7[) . we can write the zeta functions Cd^''^{s) and C)j c''('5) as 



N 



N 



i=-l 



i=-l 



where 



^A'^lc.-M = lim 



sin(7rs) 



m— >0 TT 

sin(7rs) 



vz 
a 



40F r(s + i)' 







2 \ 










— j 








. a . 




L 2 1 + z2 ; 



for i>l, 



sin(7rs) 
nm 

m— >0 TT 

sin(7rs) 



a 



a 



ztizfD'Siz)) 



dz = 



_^2s^-2s-. 



■zt(z)3M',(t(z))' 



dz = 



2s -2s— i 

—a ;^ 



r (,s + 
r(s)r(^ 



and 



pint.i./ ^ sin(7rs) 



m— >0 TT 

sin(7rs) 



z^z 
a 

vz 
a 



m ] — I log/i.(i^z) - log- 



/2^(l + z2)l/4 



AT 

E 



A(i(^)) 



dz 



X ^ [ log [c^(:^z) + vzl^{vz)\ - log 2^ — 



dz. 



Now notice that 5^(0 and can be expanded as 



hD{l) - 



l{l + D- 2)(2l + D- 2){l + D-Ay. 
{D-3)\{l + iy. 



D-2 



= J2 ( ' + 

3=0 



D-2 



boil) 



{21 + D - 2){l + D - '3)1 

^ (D~2jn\ 



D-2 



D-2 



where xd j = and yu- j = if j and D have opposite parity, and xo;o 7^ if and only if _D = 4. Let 



(A8) 



k=0 



be the Hurwitz zeta function. From (jA2p and (jASj), we find that in the case of perfectly conducting boundary 
condition, the zeta function (jAll) can be written as 



int,i+ 
D,N 



D-2 



D-2 00 
j=0 1=1 

D-2 00 ^ . , , , D-2 00 . ^ , , D-2 

,int,/+ — — , , V — ^, ,„^„int,/ + 



N D-2 

2 

i=-l 1=1 



(A9) 



E E + E (■^) + E ^'(^)^ 



i=-l 
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-,2s 



D-2 



4v^ 



{xD;j +yD;j) 



r{s + i) 



CH(2s-j-l. 



D^ 



t2s 



D-2 



{xD-j - VD-j) Ch (2s - j; y) 



T s 



-2k ^ 



i=0 

D-2 N i 

r (.s - 



_D-2 AT 



E 2^^;. E (2- + ^ - f ) E (^) 



J=0 



i+2k \ 

2 J I nint 



3 = 1 i=l 



k=0 



r(s)r(2 



where 



+ a 



1=1 

2s sin(7rs) 



2s / „-2s 



2.^ 



2 "^^l log/^^(;) (i^u(02)-log 



y2^^i;^(l + z2)l/4 ^ ^ 



AT 

E 



1=1 



- log = > ^-777- dz. 

Here I'd {I) = I + • '^'^^ result for infinitely permeable boundary condition can be obtained in a similar way. By 
taking N > D, one can guarantee that i?5^*(s) does not contain any pole on the half-plane Re s > —1. 



2. Zeta function for Casimir energy outside a spherical shell 

For Casimir energy outside a spherical shell, consider the zeta function: 

C-*(5) = C^-e*(5) + Ctm(5), 
where for perfectly conducting boundary condition, 

CXJ 

1=1 

and for infinitely permeable boundary condition, 



(AlO) 



Here 



1=1 



.^■^t,vi N v sin(7rs) 

m— >■() TT 

.„yA,vi N ,. sin(7rs) 
Cd (s) = hm 



cxt,i+- 
D 



■{s). 



^Jog[z-''K,{vz)yz, 



vz 
a 



2 .\ d 



— logjz-" [-cX,(z.z) - lyzKliiyz)]} 



dz. 



for — < Res < 1. For the Bessel function K,j{h'z) and its derivative, the Debye uniform asymptotic expansions read 
as [33, [111 : For fixed z, as z/ — > 00, 



K^iiyz) 



Klivz) 



2,(1 + ,2)1/4^^ 



E(-i)* 



■7. * ^ IJ^ 



(All) 



fe=0 
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The function Uk{t),Vk{t),ri{z),t{z) are defined in the previous section. Eq. (jAllI) implies that 

y 2 ^ 

Since 



fc=i 



log 



fc=i 



log 



^ctUk-l{t) + Vk{t) 



k=l 



(A12) 



the same approach as in the previous section shows that 

N 



N 



/"CXt,I^/ \ \ ^ *CXt.l// \ , T-)CXt,l^/ \ >cxt,i>/ \ \ ^ ^cxt,z^/ \ . T->CXt,I^ / \ 



where 



and 



\z ^int,z^ , 





= lim 

m— >0 


sin(7rs) 


/:( 


TT 






= lim 


sin(7rs) 


/:( 


TT 


d 

dz 



vz 
a 

vz 
a 



d 

dz 



\ogKt,{h'z) - log 



2i.(f +z2)l/4 



N 



V2 



E(-i) 



dz. 



1=1 



From these, we find that for perfectly conducting boundary condition, 



D-2 



D-2 AT 



r(s + i) 



J=0 



E E(-i)^c^ (2^ + * - r, f ) E ^-^^ wL r^A^ 

j=0 i=l k=0 ^ ^^'^ V 2 j 



(A13) 



E 2^^^. E(-i)^C^ (2^ + * - f ) E (^) r,,,r 

j = l i=l k=0 ^ ^^'^ V 2 i 



where 



2^sin(7rs) 



-2s / „-2s 



1=1 

TT ^ — ' 

1=1 



— I logii:^^(,) {vd{1)z) - log 



N 



72^(1 +z2)l/4 



E(-i) 



\dz 



-2s 



D-2 
2 



log 



y^^;;^g-.,,(Or,(.)(i^ ^2)1/4 N M,^o_2{t{z)) 



^/2 



E(-ir 



Mz. 



Here = / H — as before. It is easy to write down the corresponding expressions for infinitely permeable 

boundary conditions. 
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3. Zeta function for the interacting part of the Casimir energy in an annular region 



As defined in Section IIVI the zeta function tliat corresponds to the interacting term of the Casimir energy inside 
an annular region bounded by spheres of radius a and b is given by 



where for perfectly conducting boundary conditions, 



(A14) 



inter, /+ - 
D 



and for infinitely permeable boundary conditions, 



inter, /+- 



1=1 



inter, /+- 



'is). 



1=1 



1=1 



Here 



>intcr,z^/ \ i- 



sin(7rs) 



winter, z^/ \ 1- 



m— >0 TT 

sin(7rs) 



m^O TT 



(z2 -m2)-^-^log<j 1 



Iiy{az)K^{bz) 



dz. 



I^{bz)K^{az) 
[cl^{az) + azll{az)] [cK^{bz) + bzKl{bz)] 
[cl^{bz) + bzll{bz)] [cK^{az) + azKl{az)] 



dz. 



It,{uaz)Ky{vbz) / n t ti. \ t \\\ ( a(l + Vr+"Pz2-) 

- e-xjp {—2v{ri(bz) — ri(az))) — ' 



(A15) 



The Debye uniform asymptotic expansions (jA4p and (|All|) for Iu{vz) and Ki,{vz) show that as — oo. 



Iv{vbz)Ku{vaz) 



The function 



6(1 + VTT^) 



exp I 



(A16) 



a(l + Vl + 62z2) 
6(1 + Vl + a2z2) 

is < 1 for all z. Due to the power 2t^ on this function, we find that by setting m — Q directly, the first function 
on the right hand side of (|A15I) is analytic for all Res < -. It can be analytically continued to the whole plane 
using standard tricks. The possible poles of the integral appear at s = 1, 2, . . . are canceled by the zeros of sin(7rs) at 
these points. This shows that Cd'^^'^ i^) analytic function for all s. Moreover, it is obvious that as 6 — )• oo, the 

expression on the left hand side of (jA16|) vanishes. Therefore, 

hm CT'^'^is) = 0. 

o— >-C>0 

Similar argument shows that CrT'^ '^^ analytic fmiction for all s, and 

\ii^Cn.n^)-^- 

0— J-oo 



For the value of the zeta function C'"*°'^(s) at s = — -, one can immediately set s 
one finds that for perfectly conducting boundary conditions, 

Iu(i){az)K^(i){bz) 



Upon integration by parts. 



c 



inter 



{i){bz)K^^^i){az) 



1=1 ''^ 



logU 



^I^^i){az) + azllf^i^{az) 




n^K,^,^{bz) + bzKl^i^{bz)_ 




I.f^i^ibz) + bzF^^^ibz) 




^j^K^(^i){az) +azKl^^^{az) 





(A17) 



dz. 
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It is straightforward to write down the corresponding expression for infinitely permeable boundary conditions. 
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